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ABSTRACT:  A  computational  tool  was  developed  to  visualize  radar  reflection  coefficients  for  a  lossy,  two¬ 
layered  half-space  soil  geometry  as  a  function  of  incidence  angle.  Allowable  input  parameters  include  radar  fre¬ 
quency,  top  soil  layer  thickness,  and  the  complex  dielectric  properties  of  both  soil  layers.  The  tool  is  a  Microsoft 
Excel  program  that  can  operate  on  any  Windows-based  personal  computer. 
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Introduction 


A  computational  tool  is  needed  to  help  engineers  and  scientists  better 
interpret  reflectance  measurements  from  low-frequency  radars  that  illuminate 
soils.  Of  particular  interest  is  the  use  of  such  radars  to  make  nondestructive 
measurements  of  near-surface  soil  moisture  values.  Commercial  devices  are 
available  for  measuring  soil  moisture  values  that  rely  on  either  the  change  in 
response  of  a  time-domain  reflectometer  when  inserted  into  soil  (e.g.,  Dynamax 
2004;  Campbell  Scientific  2004)  or  the  change  in  shunt  capacitance  of  an  open- 
ended  coaxial  line  operating  in  the  frequency  domain  (Caldecott  et  al.  1985).  In 
both  cases,  what  is  being  measured  is  an  apparent  dielectric  constant,  which  has 
been  related  to  soil  moisture  values  through  a  calibration.  Furthermore,  these 
technologies  require  probes  being  planted  in  the  soil  or  a  person  carrying  a  probe 
to  numerous  point  measurement  locations  within  a  test  area. 

Boom-mounted  high-frequency  backscatter  radars  have  been  used  with  some 
limited  success  to  measure  dielectric  properties  of  surface  materials  and,  in  some 
cases,  to  predict  soil  moisture  values.  Ulaby’s  work  at  the  University  of  Kansas 
in  the  early  1970s  provided  an  optimistic  outlook  on  such  technology  (Ulaby  et 
al.  1974).  A  recommendation  to  use  boom-mounted  radars  operating  at  4-5  GHz 
in  a  horizontal  polarization,  horizontal  receive  mode,  and  pointing  at  the  ground 
with  an  incidence  angle  of  7-17  deg  was  quoted  in  the  literature  as  late  as  1989 
(Waite  et  al.  1984;  Schmugge  1989).  However,  Ulaby  and  his  colleagues  were 
still  not  convinced  as  late  as  1986  of  what  constituted  an  optimum  measurement 
system  design  (Dobson  and  Ulaby  1986).  Another  high-frequency  reflectometer 
was  investigated  in  the  1980s  with  applications  to  ice  that  met  with  mixed 
success  (Arcone  and  Larson  1988). 

Lower  frequency  radars  are  faced  with  the  problem  of  too  much  penetration 
of  soil  layers  under  certain  conditions  that  results  in  reflectance  values  that 
oscillate  as  a  function  of  soil  layer  thickness.  One  source  identified  a  technique 
for  using  low-frequency  radars  to  determine  layer  thicknesses  and  possibly  layer 
dielectric  properties  (Lundien  1972).  Because  there  is  strong  evidence  that 
dielectric  properties  and  soil  moisture  levels  are  correlated  (Topp  et  al.  1980; 
Curtis  2001),  one  could  theoretically  measure  dielectric  properties  and  obtain  soil 
moisture  values. 

The  objective  of  this  project  was  to  develop  a  computational  tool  that  would 
allow  one  to  easily  visualize  radar  reflection  coefficients  for  a  two-layered  soil  as 
a  function  of  incidence  angle.  Among  the  desired  input  parameters  were  radar 
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frequency,  top-layer  thickness,  and  the  complex  dielectric  properties  of  both  soil 
layers.  Such  a  tool  could  be  used  to  conduct  parametric  studies  that  show  the 
effect  of  top-layer  thickness  and  soil  properties  on  radar  reflectance  values  and 
Brewster  angles. 

Although  previous  research  has  pointed  out  some  of  the  shortcomings  of  both 
low-  and  high-frequency  radars  to  accurately  measure  soil  moisture  values,  it  is 
felt  that  there  is  still  room  for  exploring  the  phenomenology  of  the  radar  waves 
interacting  with  the  soil  layers.  The  hope  is  that  a  type  of  measurement  may  be 
discovered  that  will  uniquely  relate  reflectance  to  soil  dielectric  properties  and, 
thence,  to  soil  moisture  values.  The  computational  tool  described  in  this  report  is 
one  means  of  achieving  that  goal. 
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2  An  Approach  to  Modeling 
Layered  Lossy  Media 


Governing  Equations 

Textbooks  are  replete  with  mathematical  models  for  electromagnetic  wave 
amplitude  reflection  and  transmission  coefficients  (Stratton  1941;  Jackson  1975). 
They  account  for  complex  dielectric  properties  (although  they  usually  are  applied 
to  problems  with  real  dielectric  permittivities)  and  both  horizontally  and 
vertically  polarized  electric  fields.  They  also  cover  oblique  angles  of  incidence. 
However,  they  are  usually  limited  to  half-space  geometries.  In  other  words,  there 
is  no  allowance  made  for  multiple  layers  in  the  nontransparent  material.  If  a 
solution  is  offered  for  two  layers,  it  is  normally  done  for  the  special  condition  of 
normal  incidence. 


Half-space  horizontal  polarization 


Appendix  A  contains  the  derivation  of  the  half-space  amplitude  reflection 
coefficient  and  transmission  coefficient  relationships.  For  a  plane  electric 
wave  whose  propagation  vector  is  perpendicular  to  the  plane  of  incidence  (see 
Figure  Al),  traveling  from  material  i  to  material  j  and  striking  the  interface  with 
an  incident  angle  of  G,  the  reflection  coefficient  (the  ratio  of  reflected  amplitude 
to  incident  amplitude)  is: 


where 


ru 


£j 


o;n2  n 

CUb  C7  1 

cos#  +  —  1 

M'S/  sin2  & 

Mj  V  M 

£  <J 

— +  /■ — —  =  £•'  +  is” 

r  r*  +  •* 


(1) 


(2) 
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is  what  most  authors  agree  is  called  the  complex  relative  dielectric  constant  of 
material  j.  ^ )  is  often  called  the  dielectric  permittivity,  oj  is  the  electrical 
conductivity  of  the  material  (mho/m),  to  is  the  radial  frequency  of  the  wave,  and 
£o  is  the  permittivity  of  free  space  (8.85  x  10"12  farads/m).  p.j  is  the  magnetic 
permeability  of  the  material. 


Half-space  vertical  polarization 

When  the  incident  electric  wave  propagation  vector  lies  within  the  plane  of  "s 
incidence,  a  different  coefficient  can  be  calculated: 


cos#  -  I 

sj\ 

Mi£i 

„  £■ 
cos#  + — J 

£j\ 

M,e, 

At  normal  incidence,  with  the  wave  propagating  from  air  onto  some  nonmagnetic 
material,  both  coefficients  are  identical  and  reduce  to  the  familiar  form: 


r 


i-VJ 

1  +  yfl 


(4) 


The  vertical  polarization  equation  produces  a  negative  sign,  which  is  nothing 
more  than  an  artifact  of  how  the  electric  field  vectors  were  assumed  to  be 
pointing  in  the  derivation. 


The  three-layer  problem 

One  text  presents  a  solution  for  electric  waves  obliquely  incident  upon  a  two- 
layer  medium  (Reitz  et  al.  1979),  which  has  been  condensed  in  Appendix  B.  For 
a  plane  wave  traveling  from  material  1  (in  this  case,  air)  through  a  finite 
thickness  layer  of  material  2  and  then  into  a  semi-infinite  material  3,  the  net 
amplitude  reflection  coefficient  from  the  first  interface  is: 


r\2  +  r23g  ^ 
1  +  ri2r22e  ^ 


(5) 


where  r,2  and  r23  are  the  half-space  reflection  coefficients  at  the  first  and  second 
interface,  respectively,  and  come  from  Equations  2  and  3  for  each  polarization 
condition,  p  is  a  term  calculated  from  Equations  B1  through  B3  that  accounts  for 
wave  attenuation  through  material  2  as  well  as  a  shift  in  phase  and  is  a  function 
of  incidence  angle,  layer  2  thickness,  and  the  complex  dielectric  properties  of  the 
second  layer. 
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Obviously,  this  net  amplitude  reflection  coefficient  can  be  a  complex 
number,  indicative  of  phase  shifts  and  attenuation  within  the  top  layer  of 
material.  All  of  the  charts  that  will  be  shown  in  this  report  will  be  plots  of  the 
reflectance  of  the  waves,  the  average  energy  flux  per  unit  area.  Reflectance  will 
then  be  calculated  from  Equation  5  by  the  formula: 

R  =  rr*  (6) 


Solution  and  Display  Tool 

An  Excel  file  was  created,  a  copy  of  which  can  be  obtained  from  the  author, 
that  utilizes  these  equations  and  other  relevant  equations  from  Appendixes  A  and 
B  to  calculate  the  magnitude  and  phase  of  the  various  reflection  and  transmission 
coefficients  as  well  as  the  intensity  of  the  reflected  energy  (the  magnitude 
squared).  Once  loaded  into  a  Windows-based  computer,  the  following  procedure 
can  be  used  to  generate  results  for  specific  test  conditions: 

1.  Open  the  file  and  create  a  second  window  on  the  display  screen.  In  one 
window  of  the  display,  choose  the  “calculations”  tab.  In  the  other 
window,  choose  a  tab  for  the  type  of  chart  that  is  desired.  For  example, 
if  what  is  required  is  a  visual  display  of  the  power  reflectance  vs  angle  of 
incidence,  then  choose  the  tab  entitled  “power  reflection  (3  layers).” 

2.  On  the  upper  left  comer  of  the  “calculations”  sheet,  input  all  of  the 
relevant  operating  conditions  and  layer  properties.  For  example,  the  real 
and  imaginary  parts  of  the  complex  dielectric  constant  for  material  2 
should  be  entered  into  cells  B7  and  C7,  respectively.  Properties  for 
material  3  are  entered  into  row  9.  The  thickness  of  material  2  is  entered 
into  cell  C14,  in  units  of  centimeters.  Wave  frequency  is  entered  into 
cell  Cl 3.  Entered  data  that  effects  a  change  in  power  reflection  will 
immediately  be  seen  on  the  other  side  of  the  split  display. 

3.  A  macro  was  written  to  calculate  and  display  power  reflectance  for  any 
incidence  angle  as  a  function  of  material  2  layer  thickness.  This  chart  is 
found  at  the  tab  entitled  “power  refl  vs  layer  thickness.”  With  the 
“calculations”  sheet  active,  simply  execute  the  macro  and  view  the  result 
on  the  chart  displayed  in  the  other  window.  The  macro  has  been  written 
to  calculate  reflectances  for  layer  thickness  up  to  one  and  a  half  times  the 
wavelength  of  the  electric  field  in  material  2.  If  some  other  incidence 
angle  is  desired,  then  simply  edit  the  first  line  of  the  macro  to  set  the 
incidence  angle  value  and  run  the  macro  again. 
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3  Simulations  of  Textbook 
Problems 


To  build  confidence  in  the  computational  tool,  a  number  of  textbook 
problems  were  solved,  and  charts  of  the  results  are  presented  in  this  chapter. 
Many  of  the  published  solutions  were  for  half-space  problems,  not  multilayered 
geometries.  All  that  was  done  in  those  cases  is  that  both  layers  of  nontransparent 
material  were  given  the  same  properties.  Furthermore,  when  another  author 
presented  a  multilayer  problem,  it  most  often  was  done  for  the  special  case  of 
normal  incidence  and,  therefore,  did  not  fully  exercise  the  computational  tool. 
Nevertheless,  those  problems  present  the  only  opportunity  to  validate  the  solution 
technique  reported  herein. 

The  first  simulation  performed  using  this  computational  tool  is  that  of 
reflectance  at  an  air-glass  interface.  Results  are  shown  in  Figure  1,  which 
perfectly  match  those  found  in  numerous  electromagnetics  and  optics  texts  (i.e., 
Reitz  et  al.  1979).  Glass  is  taken  to  be  a  lossless  material  with  an  index  of 
refraction  equal  to  1.5.  Horizontal  polarization  refers  to  the  situation  in  which 
the  incident  electric  field  vector  is  perpendicular  to  the  plane  of  incidence  (for  the 
derivations  in  Appendix  A,  the  plane  of  incidence  is  the  plane  of  this  page).  For 
vertical  polarization,  the  incident  electric  field  lies  within  the  plane  of  incidence. 
The  Brewster  angle  (Equation  A50)  is  56.3  deg. 

Reflectance  charts  for  two  relatively  low-loss  natural  materials  were  found  in 
the  literature.  One,  shown  in  Figure  2,  is  for  fresh  water,  while  the  other,  shown 
in  Figure  3,  is  for  soil  (Stratton  1941).  At  100  MHz,  the  imaginary  part  of  the 
complex  dielectric  constant  for  fresh  water  equates  to  22.2  mS/m  of  effective 
conductivity.  On  the  other  hand,  the  conductivity  for  the  soil  simulation  is  1 .0 
mS/m. 

Figure  4  represents  a  simulation  of  sea  water  (reported  in  Long  1983)  at  a 
frequency  of  300  MHz.  The  effective  conductivity  of  the  sea  water  was  taken  to 
be  4.3  S/m. 

Another  interesting  simulation  taken  from  a  textbook  is  that  for  reflectance  at 
an  air-metal  interface  (Reitz  et  al.  1979).  Although  the  equations  derived  in  the 
Appendixes  are  for  dielectric  media,  this  computational  tool  was  still  able  to 
accurately  reproduce  the  published  results.  What  is  surprising  is  that  the  required 
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Three-Layer  Model 
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Figure  1 .  Reflectance  at  an  air-glass  interface 
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reflection  coefficient  amplitude  squared 
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Figure  2.  Reflectance  from  fresh  water 
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Figure  3.  Reflectance  from  a  low-loss  soil 
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Figure  4.  Reflectance  from  sea  water 
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complex  dielectric  constant  values  appear  to  be  meaningless.  What  is  meant  by 
this  is  the  following.  For  the  case  of  nonmagnetic  materials,  one  can  write  that 
the  complex  index  of  refraction  is  related  to  the  complex  dielectric  constant 
(combining  the  notations  of  Equations  All  and  A40)by 


N  =  n  +  ik  =  -yjs'  +  is” 
which  then  results  in: 

s'  =  n2  -  k 2  and  s”  =  Ink 

In  Figure  5,  the  complex  index  of  refraction  for  nickel  was  taken  to  be  2+i3. 

That  value  results  in  a  complex  dielectric  constant  of  -5+il2.  One  is  not 
accustomed  to  seeing  negative  values  of  dielectric  permittivity  when  working 
with  soils  and  other  dielectrics.  However,  negative  values  are  possible  for 
metallic  conductors  (Klein  and  Furtak  1986). 

While  this  computational  tool  is  intended  to  be  used  for  electromagnetic 
waves  traveling  from  air  to  an  air-soil  interface,  another  check  of  the  efficacy  of 
the  models  is  a  simulation  of  the  internal  reflection  problem;  i.e.,  fiber  optics.  In 
this  case,  the  conditions  are  that  the  wave  is  traveling  from  within  glass  to  a 
glass-air  interface.  Results  of  the  simulation  are  shown  in  Figure  6  and  are  exact 
up  to  the  critical  angle  of  42  deg  (Reitz  et  al.  1979).  Beyond  this  angle,  there  is 
no  energy  entering  what  is  called  the  “top  layer”  of  material  (air),  and  the 
computational  tool  attempts  to  divide  by  zero  and  fails.  Results  from  a  two-layer 
simulation  are  buried  within  the  computational  tool  and  can  be  viewed  at  the  tab 
labeled  “power  reflection  (2  layers).”  Those  results  exactly  match  the  published 
results  at  all  angles. 

A  few  textbooks  problems  involving  layered  media  are  available  as  test  cases 
for  this  computational  tool.  One  is  the  problem  of  glass  coated  with  a  dielectric 
(Reitz  et  al.  1979),  the  results  of  which  are  shown  in  Figure  7.  This  is  actually 
the  classical  quarter  wavelength  problem  for  reducing  reflectance  from  optical 
lenses.  Unfortunately,  the  textbook  authors  decided  to  chart  the  reflectance 
against  coating  thickness  normalized  to  the  free  space  wavelength,  so  that  the 
minimum  does  not  show  up  at  a  quarter  of  a  wavelength  within  the  coating. 

However,  another  author  did  show  results  for  the  quarter  wavelength 
problem  plotted  on  a  scale  where  the  destructive  interference  minima  are 
revealed  at  a  quarter  of  the  coating  wavelength  (Stratton  1941).  The  application, 
whose  simulation  results  are  shown  in  Figure  8,  is  that  of  a  dielectric  sheet  over 
water.  Perhaps  the  goal  was  to  fabricate  a  fish  tank  with  the  least  reflectance  for 
better  viewing.  As  the  results  show,  minima  do  occur  at  a  quarter  wavelength 
and  at  multiples  of  half  wavelength,  thereafter.  What  is  most  interesting, 
however,  is  that  the  simulation  shows  zero  reflectance  for  normal  incidence  when 
the  coating  dielectric  constant  is  properly  chosen.  As  displayed  on  Figure  9,  zero 
reflectance  occurs  when  the  coating  dielectric  constant  is  the  geometric  mean  of 
the  constant  for  air  and  the  constant  for  water.  Figure  10  demonstrates  that  the 
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Figure  6.  Reflectance  from  a  glass-air  interface  (total  internal  reflection) 
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Figure  7.  Normal  incidence  reflectance  from  coated  glass 
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Figure  8.  Normal  incidence  reflectance  from  a  dielectnc  coating  over  water  (e  =  2) 
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Figure  9.  Normal  incidence  reflectance  from  a  dielectric  coating  over  water  (e  =  9) 
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Figure  10.  Oblique  incidence  reflectance  from  a  dielectric  coating  over  water  (e  =  9) 
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quarter  wavelength  effect  is  still  pronounced  at  an  oblique  angle  of  incidence  but 
that  the  reflectance  does  not  go  to  zero  at  a  quarter  wavelength. 


One  final  example  is  that  of  a  sheet  of  water  with  air  on  either  side  (Stratton 
1941).  Figure  1 1  shows  the  results  of  this  simulation  using  the  computational 
tool.  It  demonstrates  two  principles.  One  is  that  reflectance  minima  will  occur  at 
resonant  conditions  set  up  at  every  half  wavelength  of  thickness.  The  second 
principle  is  that  since  the  properties  of  the  water  include  a  loss  factor,  the  peak 
reflectance  of  the  sheet  decreases  with  increasing  thickness  because  of 
attenuation  of  the  wave  energy. 
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Figure  1 1 .  Normal  incidence  reflectance  from  a  sheet  of  water 
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4  Three  Studies  on  the 
Effects  of  Material 
Properties 


Having  established  the  worthiness  of  this  computational  tool  through 
simulations  of  classical  problems  reported  in  the  previous  chapter,  it  is  now  time 
to  demonstrate  the  effects  of  material  properties  on  problems  for  which  this  tool 
was  created;  namely,  layered  soils.  The  following  pages  contain  simulations  that 
attempt  to  explore  three  different  concerns: 

1.  Are  half-space  models  “good  enough”  for  realistic  terrain? 

2.  What  is  the  effect  of  a  growing  layer  of  moist  soil? 

3.  How  is  the  response  of  a  high-conductivity  soil  different  from  that  of  a 

low-conductivity  soil? 

Dielectric  properties  for  these  soils  were  taken  from  data  collected  by  the 
author  over  a  period  of  many  years.  They  are  meant  only  to  be  typical  values. 
Soil  dielectric  properties  are  strongly  related  to  soil  moisture  and,  to  a  lesser 
degree,  soil  temperature  and  bulk  density  (Curtis  2001).  Attenuation  in  the  soil, 
represented  by  electrical  conductivity  values,  is  clearly  a  function  of  soil 
chemistry  (Curtis  2004)  and  is  virtually  unpredictable.  Unfortunately,  if  one 
wants  an  accurate  picture  of  the  soil  electrical  properties  at  any  given  site,  those 
soils  need  to  be  sampled  and  tested  in  a  laboratory. 


Half-Space  Model  vs  Layered  Soil  Model 

The  computational  tool  was  used  to  simulate  a  layer  of  moist  soil  having 
moderate  conductivity  on  top  of  a  relatively  dry  soil.  The  dry  soil  was  assigned  a 
complex  dielectric  constant  of  3  +  i0.2,  which,  at  a  frequency  of  100  MHz, 
translates  to  an  effective  conductivity  of  1.1  mS/m.  The  moist  soil’s  dielectric 
constant  was  chosen  to  be  10  +  i2,  which  results  in  a  conductivity  of  1 1.1  mS/m. 

One  question  that  begs  for  an  answer  is  “Does  soil  layer  thickness  make  a 
difference  at  a  typical  ground-penetrating  radar  frequency?”  Figure  12  is  a 
simulation  of  a  half-space  of  moist  soil,  while  Figure  13  shows  results  for  a 
5-cm-thick  layer  of  moist  soil  over  dry  soil.  Clearly,  reflectance  values  are 


Chapter  4  Three  Studies  on  the  Effects  of  Material  Properties 


19 


Three-Layer  Model 


epsilon1 

=  ( 

i 

,  0 

mu1 

=  ( 

i 

J  0 

epsilon2 

=  ( 

10 

2 

mu2 

=  ( 

i 

,  0 

epsilon3 

=  ( 

10 

,  2 

mu3 

=  ( 

1 

,  o 

frequency  (Hz) 

=  1.00E+08 

2nd  layer  thickness  (cm)  =  500.00 
2nd  layerwavelength  (cm)  =  9.44E+01 


Figure  12.  Reflectance  from  a  half-space  of  moist,  moderately  conductive  soil 
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Figure  13.  Reflectance  from  a  thin,  moist,  moderately  conductive  soil  layer  over  dry  soil 
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different  for  the  two  cases,  as  is  the  value  of  the  Brewster  angle.  The  answer, 
then,  is  “Yes!”  There  must  be  significant  penetration  of  the  thin  layer  to  produce 
such  differences. 

Fortunately,  the  computational  tool  has  the  means  to  determine  how  thick  a 
layer  must  be  before  it  appears  to  behave  like  a  half-space.  Figures  14  through 
16  show  reflectance  values  as  a  function  of  layer  thickness  at  three  different 
incident  angles.  Clearly,  for  the  properties  chosen,  the  moist  layer  would  have  to 
be  very  thick,  probably  on  the  order  of  three  wavelengths,  before  it  looked  like  a 
half-space.  Note,  also,  the  divergence  of  the  polarized  reflectances  at  higher 
values  of  incidence  angles.  Horizontally  polarized  waves  become  more 
reflective,  while  vertically  polarized  waves  become  less  reflective. 


A  Growing  Layer  of  Moist  Soil  Over  Dry  Soil 

Another  question  that  needs  to  be  resolved  through  research  is  “If  soil  layer 
thickness  makes  a  difference  on  reflectance  values,  is  there  another  measurement 
parameter  that  could  be  related  to  only  the  material  properties?”  The  rationale 
for  this  question  lies  with  the  fact  that  one  may  want  to  use  ground-penetrating 
radar  technology  to  measure  soil  properties  such  as  moisture  content.  Moisture 
content  is  clearly  a  strong  function  of  the  soil’s  dielectric  properties  (Curtis 
2001),  but  for  the  simulations  reported  above,  dielectric  properties  were  not  the 
only  parameters  that  controlled  the  response  to  the  signal.  Layer  thickness  made 
a  huge  difference. 

As  a  result  of  these  observations,  the  computational  tool  was  exercised  to 
simulate  a  growing  layer  of  moist  soil  over  dry  soil.  Some  of  the  results  of  these 
simulations  are  shown  in  Figures  17  to  21,  which  show  reflectance  vs  incidence 
angle  for  top-layer  thicknesses  of  0,  5,  10,  15,  and  20  cm.  It  is  obvious  that 
measurement  results  would  be  different  for  each  increment  of  top-layer  thickness. 
It  is  not  obvious,  however,  that  there  is  any  one  measurement  parameter  that 
remains  fixed  for  all  layer  thicknesses.  That  would  have  to  be  the  case  in  order 
for  top-layer  dielectric  properties  to  be  measurable  for  any  layer  thickness.  (As 
an  aside,  one  can  see  the  change  in  the  reflectance  curves  in  a  cartoon-like 
manner  by  displaying  the  tab  entitled  “power  reflectance  (3  layers)”  while 
executing  the  tool’s  macro.) 

The  most  obvious  measurement  parameters  that  could  possibly  be  used  to 
infer  dielectric  properties  would  be  normal  incidence  reflectance  and  the 
Brewster  angle.  Previous  results  (Figure  14)  have  already  shown  that  normal 
incidence  reflectance  oscillates  with  layer  thickness,  and  the  cartoon  simulation 
noted  in  the  previous  paragraph  shows  that  the  Brewster  angle  also  moves  around 
with  changing  layer  thickness.  As  final  proof  of  the  latter  observation, 
simulations  were  performed  for  1-cm  increments  of  thickness  all  the  way  to  a  top 
layer  of  50  cm,  which  is  a  little  over  a  half  wavelength  of  material.  Normal 
incidence  reflectance  values  and  Brewster  angle  values  were  hand  tabulated 
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Specular  Reflectance  from  Layered  Materials 
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Figure  14.  Layered  soil  reflectance  vs  top-layer  thickness  at  normal  incidence 
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Figure  15.  Layered  soil  reflectance  vs  top-layer  thickness  at  30  deg  incidence 
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Figure  16.  Layered  soil  reflectance  vs  top-layer  thickness  at  60  deg  incidence 
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Figure  17.  Reflectance  from  a  dry  soil  half-space 
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Figure  18.  Reflectance  from  a  5-cm  moist  layer  over  a  dry  soil 
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Figure  19.  Reflectance  from  a  10-cm  moist  layer  over  a  dry  soil 
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Figure  20.  Reflectance  from  a  15-cm  moist  layer  over  a  dry  soil 
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Figure  21 .  Reflectance  from  a  20-cm  moist  layer  over  a  dry  soil 
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and  plotted  on  Figure  22.  Clearly,  Brewster  angle  tracks  layer  thickness  in  the 
same  way  that  normal  incidence  reflectance  does. 

The  final  conclusion  of  this  brief  study  is  that  as  long  as  the  layered  structure 
of  soil  has  an  influence  on  the  reflectance  measurements  that  could  be  made  from 
the  soil,  there  does  not  appear  to  be  a  measurement  that  uniquely  relates  to 
dielectric  properties,  and  thence  to  soil  moisture.  The  only  hope  would  be  for 
higher  frequency  systems  where  the  skin  depth  is  much  smaller  or  for 
measurements  being  made  on  highly  conductive  soils,  which  is  the  subject  of  the 
next  section. 


Low-Conductivity  vs  High-Conductivity  Soils 

The  power  of  this  computational  tool  is  that  it  can  be  used  to  vary  one 
parameter  at  a  time  to  isolate  the  effects  of  that  parameter  on  reflectance  from 
soils.  The  first  two  sections  in  this  chapter  dealt  with  layer  thickness  as  the 
variable.  This  section  will  look  at  electrical  conductivity  of  the  top  soil  layer  as 
the  primary  variable  of  interest.  What  happens  to  reflectance  measurements  if 
the  conductivity  of  the  top  layer  changes  by  an  order  of  magnitude? 

Figures  23  through  25  are  simulation  results  for  a  moist  soil  layer  over  a  dry 
soil  layer.  Conductivity  in  the  moist  layer  was  taken  to  be  5  mS/m.  Simulations 
for  normal  incidence,  as  well  as  30  and  60  deg  of  incidence  angle  reveal  a 
behavior  similar  to  that  of  the  previous  section.  Layer  thickness  is  a  dominant 
controlling  factor,  and  polarized  signal  reflectances  diverge  for  nonzero 
incidence  angles. 

When  the  conductivity  of  the  top  soil  layer  was  taken  to  be  an  order  of 
magnitude  larger  (50  mS/m),  results  changed  dramatically.  As  seen  in  Figures  26 
through  28,  top-layer  thickness  becomes  much  less  of  a  factor  (the  skin  depth  is 
much  smaller).  In  fact,  reflectance  values  settle  into  half-space  results  for  less 
than  a  wavelength  of  top-layer  thickness.  The  second  observation  to  be  noted  is 
that  the  half-space  reflectance  appears  to  be  about  20-percent  higher  for  the  high- 
conductivity  soil  than  it  is  for  the  low-conductivity  soil. 
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normal  incidence  reflectance 


A  Growing  Wet  Top  Layer 


frequency  (Hz)  = 

top  layerwavelength  (cm)  = 

top  layer  dielectric  constant  =  ( 

bottom  layer  dielectric  constant  ~  ( 


1.00E+08 

9.44E+01 

10  >  2 
3  ,  0.2 


) 

) 


Figure  22.  Impact  of  a  growing  moist  soil  layer  on  normal  incidence  reflectance  and  Brewster’s 
angle 


Chapter  4  Three  Studies  on  the  Effects  of  Material  Properties 


brewster  angle  (deg) 


Specular  Reflectance  from  Layered  Materials 
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Figure  24.  Reflectance  from  a  low-loss  moist  layer  at  an  incidence  angle  of  30  deg 
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Specular  Reflectance  from  Layered  Materials 
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Figure  26.  Reflectance  from  a  high-loss  moist  layer  at  normal  incidence 
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Specular  Reflectance  from  Layered  Materials 
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Figure  27.  Reflectance  from  a  high-loss  moist  layer  at  an  incidence  angle  of  30  deg 
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Specular  Reflectance  from  Layered  Materials 

incidence  angle  (deg)  eo 

frequency  (Hz)  =  I.ooe+08 

top  layerwavelength  (cm)  =  8.76E+01 

top  layer  dielectric  constant  =  (  10  »  9  ) 

bottom  layer  dielectric  constant  =  (  3  ,  0.2  ) 


Figure  28.  Reflectance  from  a  high-loss  moist  layer  at  an  incidence  angle  of  60  deg 
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5  Summary  and  Conclusions 


The  objective  of  this  project  was  to  develop  a  computational  tool  that  would 
allow  one  to  easily  visualize  radar  reflection  coefficients  for  a  two-layered  soil  as 
a  function  of  incidence  angle.  Input  parameters  should  include  radar  frequency, 
top-layer  thickness,  and  the  complex  dielectric  properties  of  both  soil  layers. 

Such  a  tool  can  be  used  to  conduct  parametric  studies  that  show  the  effect  of  top 
layer  thickness  and  soil  properties  on  radar  reflectance  values  and  Brewster 
angles. 

An  Excel  file  was  created,  a  copy  of  which  can  be  obtained  from  the  author, 
that  calculates  the  magnitude  and  phase  of  various  reflection  and  transmission 
coefficients  as  well  as  the  intensity  of  reflected  energy  (the  magnitude  squared) 
for  any  given  layer  properties  and  radar  frequency.  The  validity  of  the  tool  was 
demonstrated  by  reproducing  analytical  solutions  for  layered  media  that  can  be 
found  in  numerous  electromagnetics  textbooks.  Additional  simulations  were 
performed  that  addressed  some  very  specific  situations  including  the  effect  of  a 
growing  layer  of  moist  soil  and  the  impact  of  different  soil  conductivity  values. 

One  of  the  benefits  of  this  tool  could  be  to  help  solve  the  inverse  problem  of 
deriving  soil  moisture  values  from  measured  radar  reflectance  data.  For  example, 
one  could  begin  with  the  assumption  that  a  particular  patch  of  soil  is  well- 
approximated  by  a  half-space  geometry.  Usually  one  has  normal  incidence 
reflectance  measurements  from  that  soil.  The  computational  tool  could  be 
iteratively  executed  to  find  the  complex  dielectric  constant  of  the  soil  that  best 
matches  the  measured  reflectance  value.  If  one  already  has  the  relationship 
between  volumetric  soil  moisture  content  and  the  complex  dielectric  properties  of 
that  soil,  then  a  prediction  of  the  volumetric  soil  moisture  can  be  made. 

Normally,  however,  soils  are  seldom  well-approximated  by  a  half-space 
geometry.  This  is  particularly  true  in  electromagnetic  phenomena,  where  the 
electrical  properties  of  the  medium  are  so  dependent  on  moisture  values.  The 
tool  developed  in  this  study  gives  the  user  a  means  of  simulating  more  realistic 
soil  environments.  For  example,  soils  are  never  uniformly  moist.  Soon  after 
significant  rainfall,  soils  that  exhibit  low  permeability  (such  as  clays)  may  be 
better  approximated  as  a  two-layer  half  space,  where  a  wet  top  layer  has  much 
different  dielectric  properties  than  does  the  underlying  drier  soil.  Reflectance 
values  from  the  surface  of  the  soil  will  be  a  function  of  the  amount  of  moisture 
and  the  thickness  of  that  top  layer. 
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Another  potential  solution  to  the  field  measurement  of  soil  moisture  might  be 
to  develop  a  device  for  measuring  the  Brewster  angle  (or,  at  least  the  angle  of 
minimum  reflectance)  of  the  soil  using  a  relatively  simple  bistatic  radar.  Again, 
this  computational  tool  could  be  used  to  help  solve  the  inverse  problem  of 
relating  Brewster  angle  values  to  soil  moisture  content. 
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Appendix  A 

Reflection  and  Refraction  at  a 
Plane  Interface 


Although  numerous  electromagnetics  and  optics  textbooks  present  outlines 
of  reflection/refraction  formula  derivations,  a  formal  documentation  of  the 
derivation  is  presented  in  this  appendix  that  includes  minor  details  often  assumed 
to  be  common  knowledge  by  other  authors.  Hopefully,  the  information  presented 
here  will  answer  any  questions  readers  may  have  about  assumptions  attached  to 
the  models  utilized  in  this  study.  SI  units  are  used  throughout  the  text.  The 
following  draws  heavily  upon  several  sources  (Jackson  1975;  Reitz  et  al.  1979; 
Griffiths  1981;  Schey  1973).1 


Maxwell’s  Equations 

Based  on  experimental  observations  made  during  the  early  part  of  the  19th 
century,  a  set  of  field  equations  can  be  written  that  collectively  describe  the 
interactions  among  charges,  currents,  electric  fields,  and  magnetic  fields.  Known 
as  Maxwell’s  equations,  because  of  a  final  contribution  by  the  Scottish  physicist 
James  Clerk  Maxwell  in  the  latter  part  of  the  19th  century,  they  allow  one  to 
determine  fields  from  known  charges  and  currents  as  well  as  to  predict  the  time- 
dependent  evolution  of  fields  given  a  known  set  of  initial  conditions.  In 
differential  form,  one  equation  is  referred  to  as  Gauss’s  Law: 

VD  =  p  (Al) 

which  relates  how  much  the  electric  displacement  field,  D,  diverges  (or  spreads 
out)  from  a  point  to  the  charge  density  at  that  point,  p.  In  integral  form,  this  law 
associates  the  flux  of  the  field  through  a  volume  of  material  to  the  free  charge 
enclosed  within  that  volume.  A  similar  equation  (having  no  named  source) 
accounts  for  the  observation  that  a  magnetic  field,  B,  does  not  derive  from  point 
sources: 


1  References  cited  in  this  appendix  are  located  at  the  end  of  the  main  text. 
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V-B 
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(A2) 


The  third  equation,  known  as  Faraday’s  Law,  formalizes  the  observation  that  a 
time-varying  magnetic  field  can  induce  an  electric  field,  E: 


VxE 


dB 

dt 


(A3) 


Finally,  there  is  Maxwell’s  extension  to  Ampere’s  Law,  which  accounts  for  the 
fact  that  time- varying  electric  fields  can  induce  magnetic  fields  in  materials: 


VxH 


„dD 

crE  h - 

dt 


(A4) 


where  current  density  has  been  taken  to  be  linearly  proportional  to  the  electric 
field,  and  the  constant  of  proportionality,  a,  is  the  electrical  conductivity  of 
material. 


The  following  constitutive  properties  relate  the  auxiliary  magnetic  field,  H, 
to  the  magnetic  field,  B,  and  the  electric  displacement  field,  D,  to  the  electric 
field,  E: 


H  =  —B  (A5) 

M 

and 

D  =  sE  (A6) 

where  n  is  the  magnetic  permeability,  and  e  is  the  electrical  permittivity  of  the 
material. 


The  Wave  Equation 

Assuming  that  there  is  no  free  charge  in  a  volume  of  material,  and  assuming 
that  material  is  linear  (constitutive  properties  are  all  constants),  a  wave  equation 
for  the  electric  field  can  be  derived  in  the  following  manner.  First  of  all,  apply 
the  curl  operator  to  Faraday’s  Law: 

-  -  -  -  dB 

y  x  V  x  E  —  _ V  x  — 

dt 
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Then  use  Ampere’s  Law  to  eliminate  the  magnetic  field  term: 


-  -  -  un  u  jCi 

VxVx£  =  -  ua - lie — 7- 

H  dt  *  dt2 

The  operator  on  the  left-hand  side  of  the  last  equation  can  be  replaced  by  the 
vector  identity: 

VxVx  =  VV--V2 

The  final  result  is  a  wave  equation  for  electric  fields  that  has  an  added  term 
because  of  the  conductivity  of  the  material: 

2  p  dE  d2E 

V2E  -  u<j - jus — t-  =  0  (A7) 

H  dt  dt2 


Electric  Field  and  Magnetic  Field  Solutions 

For  monochromatic  waves,  the  generalized  solution  to  the  wave  equation  is: 

E(r,t)  =  E0ei(lf-a,)  (A8) 

where  the  magnitude  of  the  wave  vector  (in  the  direction  of  propagation)  is  the 
wave  number: 

kB2£=2Ef  =  v  =  NcL  (A9) 

A  v  v  c 

N  is  the  index  of  refraction  of  the  material,  10  is  the  radial  frequency  of  the 
sinusoidal  wave,  v  is  the  wave  speed  in  the  material,  and  c  is  the  speed  of  light  in 
free  space.  Substituting  the  solution  into  the  wave  equation  gives  the  relationship 
among  the  material  properties,  the  wave  number,  and  the  index  of  refraction: 


,2  2  •  N2g)2 

k  =  fJ8(D  +1/MJC0  =  5 — 


(A10) 


or 


N2  =  fisc2  +  i^-  =  K  =  Kr+iKt  (All) 

co 

Therefore,  the  index  of  refraction  is  complex  and,  in  turn,  defines  a  complex 
dielectric  constant  whose  real  part  is: 
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Kr  =  /use1 


(A  12) 


and  whose  imaginary  part  is: 


K  =  HQ 


(A  13) 


where  fto  is  the  magnetic  permeability  of  free  space  ( 4n  X  10"7  newton/ampere2 ), 
and  So  is  the  electrical  permittivity  of  free  space  (  8.85  X  10'1"  coulomb2/newton- 
meter2 ).  For  the  vast  majority  of  materials,  the  magnetic  permeability  is  a 
number  whose  value  is  very,  very  close  to  that  of  free  space.  Only  a  few 
minerals  in  soils,  such  as  ilmenite  and  magnetite  (Telford  et  al.  1976)  have 
significant  values  of  permeability.  In  almost  all  cases,  then,  one  can  ignore  the 
ratio  of /i/juo  in  the  definition  of  the  complex  dielectric  constant. 

Given  the  solution  (Equation  A8)  to  the  wave  equation  for  electric  fields 
(Equation  A7),  one  can  easily  find  the  solution  for  the  corresponding  magnetic 
field  by  making  use  of  Faraday’s  Law.  First  of  all,  one  must  recognize  that 


k  ■  r  =  kxx  +  k  y  +  k2z 


(A  14) 


dE  _  B  jar-no  d(ik-r) 
0  cbc 


(A  15) 


Similar  terms  for  the  y  and  z  partial  derivatives  yield  the  useful  relationship  that 


V  =  ik, 

or  that  the  left-hand  side  of  Faraday’s  Law  becomes: 


(A  16) 


V  x  E  =  ik  x  E 


(A  17) 


Therefore,  the  magnetic  field  must  be  proportional  to  the  cross  product  of  the 
wave  vector  and  the  electric  field. 

As  for  the  right-hand  side  of  Faraday’s  Law,  one  can  assume  that  the 
magnetic  field  has  a  generalized  solution  like  Equation  A8  and,  therefore, 


(A  18) 


Substituting  Equations  A16  and  A18  into  Equation  A3  finally  yields: 
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1  -  - 

-kxE  (A  19) 

v 


where  k  is  the  unit  vector  in  the  direction  of  wave  propagation. 


Reflection  and  Refraction  Fields 

Following  the  notation  shown  in  Figure  Al,  one  may  now  write  expressions 
for  the  incident,  reflected,  and  refracted  electric  and  magnetic  fields: 


Incident  fields 


1  r 


E,  =  EI0ei(k'  f-M)  B,  =  — k,  x  Ej 

co 


(A20) 


Reflected  fields 


F  =  F  pr' 

‘-R 


i(kR  -r-ta) 


Br  ~  ^RX  Er 

CO 


(A21) 


Refracted  fields 


1  r 


Et  =  ETQe,(kr?-a*]  Br  =  —kT  X  E. 


CO 


J"  Z'  2-Jrp 


(A22) 


Figure  Al .  Incident  ( k ,),  reflected  (kR),  and  refracted  ( kT )  waves  at  a  plane 
interface 
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The  magnitudes  of  the  wave  vectors  are  just: 


,  N.o) 

k,  |  = 

kR 

=  kx  =  - and 

c 

kT 
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N2Q) 
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(A23) 


Boundary  Conditions 

Before  applying  boundary  conditions  to  the  field  equations,  consider  the 
benefits  of  the  assumption  that  at  any  instant  of  time  on  the  interface  boundary 
(z  =  0),  the  incident  wave  instantaneously  generates  reflected  and  refracted  waves 
having  the  same  frequency.  In  other  words,  all  three  waves  have  the  same  phase 
at  that  instant  of  time.  In  terms  of  the  wave  vectors: 


(kr-r)z=0  =  (kR-r)  z=0  =  ( kT-r)z=0  (A2A) 

Having  no  reason  to  believe  that  the  three  wave  vectors  are  not  all  in  the  plane  of 
incidence  (y-components  are  zero,  and  given  that  z  =  0  at  the  boundary,  then  the 
phase  constraint  means  that  the  x-component  of  each  wave  vector  is  identical  at 
the  boundary.  In  other  words, 

kj  sin  9l  =  kR  sin  0R  =  kT  sin  0T  (A25) 

Since  k,  =  &*  =  k,,  the  incidence  angle  must  then  equal  the  reflection  angle, 
which  we  can  call  6,.  Furthermore,  one  can  now  relate  the  incidence  angle  and 
the  transmission  angle  through  Snell’s  law: 


sin  _  k2  _ 
sin  6 2  kl  Nx 


(A26) 


In  order  to  derive  expressions  for  the  amplitudes  of  the  reflected  and 
refracted  waves  (for  it  is  the  wave  amplitude  that  is  detected  by  a  receiver 
instrument),  one  must  now  apply  boundary  conditions  for  the  electric  and 
magnetic  fields.  Given  that  there  are  no  preexisting  free  charges  on  the  interface, 
nor  are  there  any  preexisting  currents  on  the  interface  because  of  those  charges, 
then  in  words  one  can  state  that  the  normal  components  of  the  electric 
displacement  field,  D,  and  the  magnetic  field,  B,  must  be  continuous  across  the 
interface,  as  are  the  tangential  components  of  the  electric  field,  E,  and  the 
auxiliary  magnetic  field,  H.  In  equation  form,  the  boundary  conditions  are: 


[^1  C^/0  ER0  )  £2^To\’^  0 

(A27) 

| kj  *  EI0  +  kR  x  Ero  —  kj  x  E^q  ]•  n  —  0 

(A28) 

o 

+ 

>3 

O 

1 

> 

o 

X 

S> 

II 

o 

(A29) 

A6 


Appendix  A  Reflection  and  Refraction  at  a  Plane  Interface 


0 


(A30) 


(kj  x  Ej o  +  kR  x  Ero  )  ( kT  x  iij.0 

Mi  Mi 
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Amplitude  Reflection  and  Transmission 
(Refraction)  Coefficients 

Because  any  incident  electric  field  vector  can  be  defined  by  two  orthogonal 
vectors,  amplitude  reflection  and  refraction  coefficients  (ratios)  will  be  developed 
for  two  special  cases:  one  in  which  the  electric  field  is  normal  to  the  plane  of 
incidence  (horizontal  polarization),  and  one  in  which  the  field  lies  within  the 
plane  of  incidence  (vertical  polarization). 


Horizontal  polarization 

For  this  special  case,  the  field  amplitude  vectors  and  wave  propagation 
vectors  take  on  the  following  forms: 


EI0 

A  — *.  A  —  A 

=  El0  j  ,  Er o  =  ErJ  ,  Eto  =  ETOj 

(A31) 

kj 

=  ^  sin^/  +  kx  cos^^ 

(A3  2) 

kR 

A  A 

=  kx  sin  0X  i  -  kx  cos  6X  k 

(A33) 

kT 

A  ^ 

=  k2  sin  02  i  +  k2  cos  02  k 

(A34) 

Inserting  these  vector  definitions  into  the  first  boundary  condition,  Equation  A27, 
yields  nothing.  The  second  boundary  condition,  Equation  A28,  combined  with 
Snell’s  law,  yields  the  same  relationship  that  one  gets  from  the  third  boundary 
condition;  namely: 

En  +  Er  o  —  Et  0  =  0  (A3  5) 

The  fourth  boundary  condition,  Equation  A30,  yields: 

k  k  k 

— —EI0  cos  9X  A — -ER0  cos  0X  h — —Eto  cos  92  =  0  (A36) 

Mi  Mi  Mi 

Equations  A35  and  A36  can  be  solved  simultaneously  to  arrive  at  a  reflection 
amplitude  coefficient  and  a  refraction  amplitude  coefficient.  For  example, 
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cos#,  -  cos 02 

_ MiK 

COS  6,  +  cos  Q 

MiK 

Known  parameters  will  be  the  incidence  angle  and  the  optical  properties  of  each 
layer.  The  transmission  angle  can  be  eliminated  from  Equation  A37  by  using 
both  Snell’s  law  and  the  trigonometric  identity:  sin2  0  +  cos2  0=1. 
Furthermore,  the  wave  number  in  each  layer  can  be  replaced  by  terms  that 
include  either  the  complex  index  of  refraction  or  the  complex  dielectric  constant 
defined  by  Equation  All: 

cos  6X  -  — 

Ero  _  _ th. 

Eio  a  ,  A 

COS  +  — 

Mi 


Using  Equations  A12  and  A13  to  express  the  reflectance  amplitude  coefficient  in 
terms  of  dielectric  permittivity  and  electrical  conductivity,  one  finally  has: 


(A3  9) 


Other  authors  report  the  same  relationship  shown  in  Equation  A39,  but 
usually  express  it  in  terms  of  a  complex  index  of  refraction.  This  author’s  recent 
body  of  work  has  exclusively  dealt  with  measurements  of  the  complex  dielectric 
properties  of  soils  under  the  assumption  of  nonmagnetic  behavior.  In  that  case 
the  complex  relative  dielectric  constant  for  layer  1  is  just 


ex  =  +  =  e[  +  ie”x  (A40) 

£0  £q(0 

and  Equation  A3  9  can  be  expressed  in  the  form  inserted  into  the  Excel  file  used 
to  generate  the  reflection/refraction  results  shown  in  this  report: 
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EI0 

cos#!  +—  I 

MiS2 

-  sin2  #, 
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Combining  Equations  A3  5  and  A41  yields  the  corresponding  transmission  (or 
refraction)  amplitude  coefficient: 


2cos^ 

cos#!  +  —  i  ^2f2  -  sin2  #, 
Mi  V  M\e i 


(A42) 


Vertical  polarization 

For  the  case  of  vertical  polarization,  all  of  the  electric  field  vectors  lie  in  the 
plane  of  incidence  and  take  the  forms: 


EI0 

yy  A 

=  -  EJ0  cos  #!  i  +  EI0  sin  #,  k 

(A43) 

Ero 

=  Ero  cos  #,  i  +  Em  sin  #,  k 

(A44) 

Eto 

A  A 

=  -  Eto  cos  #2  i  +  Eto  sin  #2  k 

(A45) 

The  wave  propagation  vectors  are  identical  to  Equations  A32  through  A34. 

Applying  these  vector  definitions  to  the  first  boundary  condition  (Equation 
A27)  yields  an  equation  identical  to  that  derived  by  applying  the  fourth  boundary 
condition  (Equation  A30): 


mA\ 

(  E  i 
1+  R0 

M\£i 

F 

\  o  2 

(A46) 


In  this  case,  all  of  the  substitutions  described  above  for  the  case  of  horizontal 
polarization  have  already  been  applied  to  derive  the  transmission  amplitude 
coefficient  in  terms  of  the  reflection  amplitude  coefficient. 

The  second  boundary  condition  (Equation  A28)  provides  nothing,  while  the 
third  boundary  condition  (Equation  A29)  gives: 


cos#!  f  _  er; 

cos#2^  £/0y 


(A47) 
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Setting  Equation  A46  equal  to  Equation  A47  provides  a  formula  for  the 
reflection  amplitude  coefficient  for  vertical  polarization  in  terms  of  layer 
properties  and  the  angle  of  incidence: 


COS0.-J 
^ RO  _  1 

M  2&l 

(e  '] 

b\ 

2 

sin2  9X 

M\£2 

V.^2  J 

F 

^7  0 

COS0,  +  J 

2 

sin2  6X 

f*\S2 

J 

(A48) 


Finally,  substituting  this  result  into  Equation  A46,  one  has  a  relationship  for  the 
transmission  amplitude  coefficient  for  vertical  polarization: 


2cos#! 


^^fcos  9X  + 

i 


1- 


/fifi 


sin2  Qx 


(A49) 


Brewster’s  Angle  and  the  Critical  Angle 

For  lossless  materials,  there  are  conditions  for  which  reflectance  is  zero.  The 
condition  is  trivial  for  horizontal  polarization  problems.  Setting  the  numerator  of 
Equation  A41  to  zero  and  assuming  nonmagnetic  properties,  one  has: 


cos  <9, 


-  sin2  9X 


Squaring  both  sides  of  the  equation  and  using  one  of  the  fundamental 
trigonometric  identities,  one  has: 


In  other  words,  the  two  materials  are  identical;  the  interface  does  not  exist. 

For  vertical  polarization,  setting  the  numerator  of  Equation  A48  to  zero  and 
applying  the  same  assumptions  and  substitutions,  one  gets: 

COS  0  =  - - y- 

*2  “*1 

Inverting  the  equation  and  substituting  the  trigonometric  identity  involving  sec2 
and  tan2,  one  finally  gets: 
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The  angle  given  by  this  equation  is  called  the  Brewster  angle  and  represents  the 
incidence  angle  at  which  no  vertically  polarized  energy  can  be  reflected.  As  an 
aside,  since  electromagnetic  waves  can  be  viewed  as  a  linear  combination  of 
horizontally  and  vertically  polarized  components,  Brewster  angle  reflection  can 
be  a  tool  for  producing  linearly  polarized  energy. 

If  one  then  substitutes  Snell’s  law  into  Equation  A50,  one  finds  that  at  the 
Brewster  angle, 

cos#fl  =  sin#2 
which  can  only  be  true  if 


Therefore,  at  the  Brewster  angle,  the  angle  of  incidence  and  the  angle  of 
refraction  are  complementary. 

There  are  also  conditions  for  which  reflectance  is  total.  The  incident  angle 
for  which  those  conditions  are  met  is  called  the  critical  angle.  One  can  find  the 
critical  angle  either  by  setting  the  reflection  amplitude  formulae  (Equations  A41 
and  A48)  to  unity  or  by  setting  the  refraction  angle  in  Snell’s  law  to  90  deg.  In 
either  case,  the  critical  angle  is  found  from: 


sin<9c 


(A51) 
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Appendix  B 

Reflection  and  Refraction  from 
a  Thin  Slab 


The  following  figures  and  formulae  have  been  lifted  directly  from  a  popular 
electromagnetics  theory  textbook  (Reitz  et  al.  1979).  The  derivations  are 
extensive,  and  interested  readers  are  encouraged  to  read  the  source.  This 
appendix  serves  only  as  documentation  for  the  equations  used  to  perform  the 
multilayered  wave  propagation  simulations  reported  in  earlier  sections. 


Phase  Shifts  and  Amplitude  Attenuation 

Two  figures  are  required  to  visualize  the  problem  of  thin  slab  reflection  and 
refraction.  The  first  is  useful  in  understanding  how  to  account  for  the  phase 
differences  in  energy  that  emanates  from  a  plane  interface  because  of  multiple 
reflections  and  refractions  of  waves  inside  of  a  finite  thickness  layer  of  material. 
The  situation  is  this:  A  plane  wave  traveling  through  a  material  defined  by  an 
index  of  refraction,  Nj ,  impinges  upon  a  slab  of  material,  N2 ,  whose  thickness  is 
d.  One  can  think  of  two  parallel  rays  in  the  direction  of  the  wave  propagation 
vector  in  Ni,  one  of  which  strikes  the  slab  at  O  and  the  other  which  strikes  at  X. 
The  incidence  angle  is  0i.  The  ray  at  O  partially  reflects  and  partially  refracts 
into  N2.  The  amplitudes  of  the  reflected  and  refracted  waves  are  determined  by 
the  formulae  developed  in  Appendix  A.  The  refracted  portion  of  the  incident 
wave  travels  through  material  N2  where  it  reflects  off  and  refracts  into  material 
N3.  Again,  this  behavior  is  controlled  by  the  formulae  in  Appendix  A.  The 
reflected  portion  passes  back  through  N2  and  partially  reflects  back  into  the  slab 
and  partially  refracts  into  material  Ni.  The  wave  that  comes  out  of  N2  has  not 
only  been  attenuated  by  the  material  within  the  slab,  but  it  has  experienced  a  shift 
in  phase  relative  to  the  incident  ray  that  impinged  upon  the  slab  at  X.  The 
reduction  in  amplitude  and  the  shift  in  phase  can  both  be  accounted  for  by  an 
exponential  term,  eifi,  where: 

P  =  2  d—(p  +  iq)  (Bl) 
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The  real  part  of  ft  gives  the  phase  shift,  while  the  imaginary  part  gives  the  two- 
way  attenuation  of  the  wave  amplitude  within  the  slab.  In  terms  of  the  incident 
angle,  the  dielectric  constant  of  Ni  and  the  complex  dielectric  constant  of  the 
slab,  the  components  of  ft  are  found  from: 


Figure  B1 .  Paths  of  two  rays  incident  on  a  slab  of  finite  thickness 
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These  equations  are  the  same  as  those  that  relate  the  components  of  a  complex 
dielectric  constant  to  the  components  of  a  complex  index  of  refraction,  except 
that  K2r  is  replaced  by  (K2r  -  Kisin  0!).  In  other  words,  p  and  q  are  generalized 
versions  of  the  real  and  imaginary  parts  of  the  complex  index  of  refraction  that 
depend  on  the  angle  of  incidence. 


Multiple  Reflections 

A  second  figure  is  useful  to  visualize  the  geometry  of  multiple  reflections 
within  the  slab  that  will  lead  to  the  final  formulation  for  reflection  and 
transmission  amplitude  coefficients. 


Figure  B2.  Reflection  and  refraction  of  a  single  ray  incident  on  a  slab  of  finite 
thickness 
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Making  use  of  the  notation  that  ry  is  the  amplitude  reflection  coefficient  for  a 
wave  traveling  from  material  i  to  material  j  (either  Equation  A41  or  Equation 
A48),  and  ty  is  the  amplitude  transmission,  or  refraction  coefficient  for  a  wave 
traveling  from  material  i  to  material  j  (either  Equation  A42  or  Equation  A49), 
one  can  write  down  an  algebraic  series  that  adds  up  all  of  the  contributions  to  the 
net  amplitude  reflection  coefficient,  r. 


r  —  ^12  ^12^23^21^23^21^  ' 


(B4) 


—  rn  +  hiriJi\e 


iP 


|1  +  7*21  *23e  ^  i,r2ir23e  P  'f 


+  . 


Using  the  infinite  series  identity: 


l  +  z  +  z2  + 


1 

1  —  z  ’ 


one  can  rewrite  Equation  B4  as: 


=  r\2  + 


t  t  r  p'P 

f12f21A23c 

1  —  Y  Y 
1  r2Y23c 


_  r\2  +  ^23  (^12^21  ri2r2\ 

1  ~  V2\r23eP 

Finally,  making  use  of  the  identities: 


'12 


=  ~r- 


21 


(B5) 


(B6) 


and 


^12+^2.  =  1>  (B?) 

one  obtains  the  final  form  of  the  amplitude  reflection  coefficient  as 
Y  +  Y 

r  _  M2  (B8) 

1  +  r\2r23e  ^ 

The  same  procedure  could  be  used  to  find  expressions  for  the  amplitude 
transmission  coefficient  into  material  2,  the  amplitude  reflection  coefficient  for 
the  interface  between  material  2  and  material  3,  and  the  amplitude  transmission 
coefficient  into  material  3.  However,  the  only  coefficient  needed  for  this  study  is 
that  given  by  Equation  B8. 


B4 
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It  should  be  noted  that  these  relationships  apply  for  either  horizontal  or 
vertical  polarization.  The  two-material  reflection  coefficients  in  Equation  B8 
come  from  the  expressions  derived  in  Appendix  A  (either  Equation  A41  or 
Equation  A48),  and  the  phase  factor,  /?,  comes  from  Equations  B1  through  B3. 
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